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We show that a recently proposed [J. Fleischer and P.H. Diamond, Phys. Rev. E58, R2709 
(1998)] one-dimensional Burgers-like model for magnetohydrodynamics (MHD) is in effect identical 
to existing models for drifting lines and sedimenting lattices. We use the model to demonstrate, con- 
trary to claims in the literature, that the energy spectrum of MHD turbulence should be independent 
of mean magnetic field and that cross-correlations between the noise sources for the velocity and 
■ - - , magnetic fields cannot change the structure of the equations under renormalisation. We comment 

' on the scaling and the multiscaling properties of the stochastically forced version of the model. 
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The search for simple model equations embodying some of the features of fluid turbulence prompted Burgers to 
propose his famous nonlinear diffusion equation H 
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I. INTRODUCTION AND RESULTS 



in one space dimension, where fo is a "viscosity" and / a forcing function. The unique properties of (|^), particularly 
the Cole-Hopf transformation that linearises it, together with its obvious formal similarity to the Navier-Stokes 

■ equation, have led to many studies shedding some light on real fluid turbulence. Turbulence in plasmas is believed 
to be described by the equations of magnetohydrodynamics Q| in three dimensions (3dMHD) for the coupled evolution 

Ci I of the velocity field u and the magnetic field B. In SdMHD, the Navier-Stokes equation for an incompressible fiuid is 
^ modified by the inclusion of electromagnetic stresses: 

^ ' du ^ Vp (V X B) X B _o , 

■ — + u.Vu = + ^ + iyV^u + fu 2 

5 ■ dt p 4ttp 

O 



with V.u = 0; and Ampere's law for a conducting fluid becomes 

— U.VB = B.Vu -H^V^B + fb. (3) 
ot 

In (||) and (||), p is the fluid density, p is the pressure, p is the "magnetic viscosity" , arising from the nonzero resistivity 
of the plasma, j/ is the kinematic viscosity and fu and fb are forcing functions. Since the MHD equations are distinctly 
more complicated than the Navier-Stokes equation, a toy model, whose relation to the MHD equations is the same 
as that of the Burgers equation to the Navier-Stokes equation, should be very welcome. 

The first attempt in this direction |^ lacked cross-helicity conservation. The Id model with all the sca/ar conservation 
laws of 3dMHD is that due to Fleischer and Diamond ||]. Their equations, however, have been in the literature for 
some time now in a different context. For the case where the mean magnetic field = 0, the equations are those 
of the Erta§-Kardar model of drifting lines. For the general case bo ^ 0, they are the coupled equations, for 
the displacements along and normal to the drift direction, of the one-dimensional reduced model of Lahiri and 
Ramaswamy |^ for a crystalline lattice moving through a dissipative medium. This surprising relation arises from 
an exact correspondence between the symmetries of these equations and those of the IdMHD equations with and 
without bo — 0. The relation between the velocity and magnetic fields in Q and the displacement fields in is 
precisely the same as that between the height in the KPZ equation and the velocity in (Q). 

Reduced models are useful if they allow one to answer a question of principle. We use the IdMHD model to show 
that the form of the energy spectrum in MHD with a mean magnetic field should be the same as that without one, 
contrary to a claim by Kraichnan [p^ . We also explore the scaling and multiscaling properties of the model equations 
in the presence of stochastic forcing. For a temporally white noise source with variance oc /c"^"*"^ at small wavenumbers 
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A;, we find "sweeping" divergences for y > 3 and, within perturbative and renormalisation-group analyses, we show that 
the effective kinematic and magnetic viscosities are equal at long wavelengths. Specifically for y = 4, as is appropriate 
for a study of Id turbulence, we find within a scaling treatment that the energy spectrum E{k) ~ k~^/^ and the 
dynamic exponent z = 2/3. We show also that the multiscaling properties of the model should be precisely the same 
as that of a Burgers equation with the same driving . Lastly, we correct a claim jl^] in the literature regarding the 
role of cross-correlations between the noise sources in the velocity and magnetic-field equations. 

II. SYMMETRIES, EQUATIONS OF MOTION, AND NOISE STATISTICS 

Let us remind the reader of the "derivation" of the IdMHD equations Q emphasising some points which we feel 
are important and are missing in ref. Q. The structure of the Id model is determined by the invariances of the 
MHD equations, (i) The velocity and magnetic fields are polar and axial vectors respectively. For a model in one 
space dimension with one-component fields u{x) for the velocity and b{x) for the magnetic field this implies evolution 
equations which are invariant under x — s- — x together with u — s- —u, with h b. (ii) The MHD equations are Galileian- 
invariant, so our Id model must be unchanged under x ~ x' + uqI' , t — t', u{x,t) = u'{x',t') + uq, b{x,t) = b'(x',t') 
for any constant uq. (iii) The Id equations, like the 3d originals, must have the "conservation law" form dt{) = V.( ). 
(iv) The equations, in the absence of the diffusion terms, must conserve the "energy" 



{u^+(3b^)dx (4) 

for some constant /?. To leading orders in gradients and bilinear order in fields, the most general Id MHD model 
consistent with these requirements Q is Q 

du . , db du , , db d'^u . , 

+ Mbo-^ + Mu— + X2b— ^ Vo—^ + U, (5a) 
ot ox ox ox ox'^ 

9b , ^ d , a , d^b ^ 

m + ^^^"a^" + ("^^ = + ^^^^ 

Integrating (||) once with respect to x, one finds that u and b obey the equations of Erta§ and Kardar 0| 
for drifting lines when bo — 0, and those of Lahiri and Ramaswamy |^ for sedimenting lattices when bo ^ 0. This 
transformation of variables is the same as that which takes one from the Burgers equation to the KPZ ^ equation. 

In (|^), the constant bo is the "mean magnetic field", v and /i are the kinematic and magnetic viscosities, and 
/„ and fb are forcing functions. Ai can always be rescaled to unity but is retained for book-keeping purposes. A2 
is arbitrary, and all statements made in this paper hold independent of its value . It is straightforward to show 
that the equations, in the absence of the diffusion terms, conserve the energy- like quantity £ defined in eq.(|4|), for 
/? = A2/A1, as well as ub, the analogue in d = 1 of the cross helicity (J^ u.B) [Q, and that the model of does not 
conserve the latter. A third conserved quantity, the magnetic helicity / A.B, where A is the vector potential, has no 
analogue in d = 1. It is important to note that for bo=0 the equations have a higher symmetry than for 60 7^ 0. For 
both 60 = and bo ^ 0, they are of course invariant under x — s- — x, u — > — u, 6 — > 6, as imposed by the transformation 
properties of the fields. For bo = 0, they are in addition invariant under x —x,u —u, b — > —b. For A1A2 > 0, in 
the linear approximation at small wavenumber k, as discussed in S, the equations (0) have a wavelike response with 



frequency u) = ±a/AiA2^. These are the Id analogue of Alfven waves [||J|]. The hnearly unstable case A1A2 < ||] 
will not be considered here since the corresponding possibility is unphysical in the context of the 3dMHD equations 
(iandl). 

Since we wish to study turbulence in (|^), we shall take the forcing terms /„, /{, to be zero-mean Gaussian noise 
sources with covariance specified below. The properties of u and b under x — > —x and the reality of fuix, t) and ftix, t) 
imply that < fu{k, 0)fu{—k, t) > and < fb{k, 0)/f,(— fc, t) > are real and even in k while, crucially, < fu{k, 0)fb{~k, t) > 
(if nonzero) is odd in k and purely imaginary. In general, therefore, 



^The choice of bookkeeping coefficients in the equations as presented in Q seems unnatural to us. Galileian invariance requires 
equality of the coefficients of udxU and dx{ub), while that of hd^b can be left entirely arbitrary. 
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</n(fc,0)/„(fc',t) > = A^{\k\)5{k + k')S{t) (6a) 
<ft{k,0)h{k',t)> = At{\k\)Kk + k')S{t) (6b) 
< Uk,0)hhk,t) > =ikC{\k\)S{k + k')6{t), (6c) 



where A^, At, and C are real functions. By strict analogy with the random stirring approach in d = 3 ||12|,|13[ and Id 
Burgers turbulence |^ we choose 

Mk) = ^ (7) 

so that the noise strength has the units of the dissipation parameter, i.e. (length)'^/ (time) (fufb) and {fbfb) can 
be taken to be zero or at any rate no more singular than (fufu) without altering any of our main conclusions. A 
self-consistent treatment starting with only a nonzero (fufu) will generate (fufb) and (fbfb) as well. In that sense, 
contrary to the statements in the case where only one of the equations has a bare noise is not physically distinct 
from that where both noises are non-zero. For convenience in carrying out bare perturbation theory, however, we take 
all components of the noise covariance to be as singular as (fufu)- 

M\k\) = ^ (8) 
kCi\k\) = '-f (9) 

Note that it is 1/fc and not l/\k\ that appears on the right-hand side of (^. We now use the reduced equations (||) 
to obtain the results advertised in the Introduction. 

III. ENERGY SPECTRUM AND THE ROLE OF A MEAN MAGNETIC FIELD 

Does the presence of a mean magnetic field 6o affect the energy spectrum E{k) at wavenumber k in MHD turbulence? 
It has been argued |ic| ] that if bp 7^ then (i) the dissipation e must be proportional to the Alfven wave propagation 
time 1/cfc, where c is the wavespeed; (ii) e must depend, apart from this, only on k and E{k)] and (iii) as a result of 
energy conservation and a local-cascade picture, e must be independent of k. Taken together, these imply [ic[| that 
E{k) ~ fc~^/^, while for = the usual Kolmogorov arguments yield 

The questionable assumption in the foregoing analysis is that the wave-propagation time sets the scale for the 
dynamics. In what follows, we test that assumption in the IdMHD (^ model since the above arguments, if correct, 
should apply there as well. We show below explicitly that the presence of a nonzero wavespeed leaves the equal-time 
correlations unchanged at small k. We begin by rewriting (||) in terms of the Elsasser |14| variables 



z± = , ■ (10) 

As has already been noted in 0,^, the equations for z+ and z~ decouple completely |^ at small wavenumbers 0: 



— -6o5:,z+ + -^-^ d^z+ ^ Dd^^z+ + ... + f+- (11a) 

dz~ 
'dt 



' bod^z +\ ^^^'^ d^z ^ = Dd^^z +... + / , (lib) 



with 



/„ ± VX2fb 



(12) 



■^The decoupling is demonstrated in |^ for zero mean magnetic field, with equal magnetic and kinematic viscosities. Our 
treatment is somewhat more general. 
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where we have set Ai to unity, D = (/z + i/)/2, and the elHpsis refers to reactive (i.e., nondissipative) terms proportional 
to V — fi^ subdominant in wavenumber relative to the leading Alfven wave terms (x bo. It therefore suffices to study the 
case V = ^1 = D . In (Ila) and ( |llb| ), the wave terms cx 60 can be absorbed by opposite Galileian boosts, x ^ x — bot, 
and X X -\- hot respectively. Each equation can therefore be studied independently, at small wavenumber, as far 
as the autocorrelations of or z~ are concerned. However, the correlations of the physical fields u and b involve 
correlations (z+z^) (which are nonzero because {f^f^) 7^ 0). Some care must be taken while evaluating these, as we 
will show in detail below. Since 



7\ (z^ + z ) (13) 



and 



we see that 



b=^q^iz^-z-) (14) 



(u{k,t)u{^k,0)) = ^-^[{z+{k,t)z+{-k,t)) + {z-{k,t)z-{-k,t)) +2Re{(z+{k,t)z-{~k,0)))] (15) 
{b{k, t)b{-k, 0)> = 1^ [(z+(fc, t)z+{~k, t)) + (z-{k, t)z- i-k, t)) - 2Re((z+(fc, t)z-{~k, 0)))] (16) 



and 



{u{k,t)b{-k,0)) = l^thl[{z+{k,t)z+{-k,t)) - {z-{k,t)z-i-k,t)) -21mi{z+{k,t)z-i~k,0)))] (17) 
4VA2 

Let us define shifted coordinates 

y±^xT bot (18) 

and "comoving" fields 

C±(y±,i) = C^ixTbot,t) = z^ix,t). (19) 

Then, Fourier transforming, we get 

z=^(fc,i) = e=F*'^''°*C'^. (20) 
The autocorrelations of z+ and C+, as well as those of z~ and C~, are related by travelling waves: 

(z±(fc,i)z±(fc',t')) = eT^'=^°(*-*')(C±(fc,i)C^(fc',0), (21) 
while that between z+ and z^ has a phase factor which depends on the sum of the time arguments: 

{z+{k,t)z-{k\t')) = e-'^''°(*-*')(C+(A;,i)r(fc',t')>- (22) 

Defining 

A ^ (23) 
we see that the fields obey the Burgers equation 

dtC± + ^dyCi = DdyyCi + <P±{y, t) (24) 
with noise correlations related to those of /+ and in eq. ( |l^ ) by 

{<t>+{k,t)^+{k',t')) = (/+(A:,i)/+(fc',t')), (25) 
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{r{k,t)r{k',t')) = (/-(fc,t)/-(fc',o), (26) 

and 

{cf>+{k,t)r{k',t')) = e*'='"'(*+*')(/+(fc,i)/-(fc',0)- (27) 

Despite the appearance of nonstationary phase factors in the correlations {4)'^(t)~) and (C^C~)i ^ result of the time- 
dependent coordinate transformation, all physical correlation functions will of course be time-translation-invariant. 
The autocorrelations (C^C^) ^^id (C~C~) particularly simple and entirely independent of the wavespeed 6o- Thus 
the energy spectrum which, for Ai = 1, is 

^^^^^\ J ^[{Hk)f)+X,{m\')] (28) 

= ^^/^[(k+(fc)P> + (r(fc)P>] (29) 

= '-^J^mik)\') + {\C{k)\')] (30) 

is therefore identical to that of the randomly forced Burgers equation ( p^ ) , which does not contain the Alfven waves 
at all. This refutes the claim of since the arguments in that work would, if correct, have applied to the present 
model as well. This establishes the assertion in the Introduction that a nonzero mean magnetic field is irrelevant to a 
determination of the energy spectrum of MHD turbulence. It is worth noting that our result holds regardless of the 
nature of the forcing (deterministic, stochastic, with or without long-range correlations). 

A demonstration that equal-time correlations in ^ are independent of bo, requires showing that {z~^{k)z~ (—k)) 
does not involve the wavespeed. We establish this below for sufficiently small wavenumber. We must assume in our 
derivation that the long-wavelength behaviour of (E4) is determined by a fixed point at which Galileian-invariance 



prevents the nonlinear coupling in (24) from renormalising, and at which higher-order nonlinearities are not relevant. In 
that case, the renormalised {z^ {k)z~ {—k)) correlation function can be written in terms of the renormalised propagator 

G{k,t) of the randomly-forced Burgers equation ( p^ and the (renormalised, in principle) correlator {k,t) of 

f^{k,t) with /^(— fc, 0) with appropriate phase factors arising from the change of variables ([T^): 

{z+{k,t)z-ik',t')) S{k + k') [ dh f dt2G(fc,t-ti)iV+_(fc,ti -i2)e*''°''(*^+*')G(-A:,i'-i2) (31) 

= S+-{k,t-t')5{k + k') (32) 
Redefining t ~ ti ~* ti, t' — t2 — > ^2, we obtain the equal-time correlations 

/"OO poo 

S+-{k)= dti dt2G(/c,ii)G(-fc,t2)e*'=''°(*i+*=)iV+_(/c,t2 -ti) (33) 
Jo Jo 

Scaling tells us that 

G{k,t) = T{kH) (34) 

N+^{k,t) = k^R+^{k)TN{kH) (35) 

where z is the dynamic exponent for equations (p^, i.e., without the wave terms, F and Fat are scaling functions, and 
ik) is the zero-frequency renormalised covariance of and Defining ri = k'ti and T2 — fc^i2, we see that 

/•oo i*oo 

S+-{k) = R+.{k) dri dr2^(Tl)^(r2)^^^(rl-T2)e'(^^+^^)''°'='"^ (36) 
Jo Jq 

Since the dynamic exponent governing the decay of correlations in the Burgers equation with noise variance ~ 1/fc 
is known to be 2/3, the phase factor inside the integral ( ^6| ) for small k can be set to unity, i.e., the wavespeed 
term drops out. This completes the demonstration that equal-time correlations in MHD turbulence are independent 
of the mean magnetic fieldR. 



and 



^Even if we had a model with z > 1, say, for a noise with vanishing variance at zero wavenumber, the contribution from 
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IV. SCALING, RENORMALISATION GROUP AND SELF-CONSISTENT ANALYSES FOR ZERO 

MEAN MAGNETIC FIELD 



Some general observations will allow us to estimate the scaling exponents for this IdMHD model. First the Galileian 
invariance mentioned before implies that Ai will not renormalise. Secondly, the fact that bo can be eliminated 
completely from equation ( |lO| ) to give equation s (p^ by the transformations x —> x zL bot means that bo cannot 
renormalise. Thirdly, the invariance of equations (|24|) under the Galileian transformation C^{x, t) = (^[x + X2Uot, i)± 
Uo, means that the coupling strength A2 in equations (|2j) cannot renormalise 0. 

Thus, if we were to carry out a renormalisation-group transformation by integrating out a shell of modes Ae~^ < 
q < A, and rescaling x — > e^x, u e^^^u, b e^^''b, t e/^t), the couplings Ai and A2 would be affected only by 
the rescaling: 

Ai ^ e^(X"+^-i)Ai, A2 e^(2x.-x.+-i)A2. (37) 

We can thus rescale to keep Ai and A2 fixed, giving 

X« = Xb = 1 - 2. (38) 

The noises sources /„ and fb have variances (see equation(|^)) which diverge at small wavenumber k. Since the nonlinear 
couplings are first order in k, any renormalisation of eu,€b and e„f, will be of the form k'^W{k,uj) which vanishes for 
fc — > 0,w 7^ 0. The noise is thus unrenormalised at small k for any nonzero frequency. If we extend this to say that 
the noise strength receives no fluctuation corrections at all, then, under a renormalisation-group transformation, the 
noise strength too is affected only by rescaling. Insisting that the rescaling leave the noise strength unchanged yields 
z = 2xu = 2x6- Thus, z = 2/3 and Xu = Xb = 1/3, so that E(k) ~ While this treatment, which is the same 

as that applied to the Burgers equation with noise variance 1/fc, does not yield niultiscaling, it appears to be 
satisfactory for two point correlations. Since our equations for (|2^ ) are identical to jsj, it is reasonable to expect 
that the IdMHD equations (|^) with singular noise will show multiscaling similar to that in In particular, in this 
Id treatment, z^ and z^ obey the same equation implying that their multiscaling properties are the same. This is 
consistent with the behaviour of the Elsasser fields in a shell model of 3dMHD |l5j . 

In section (III) our decoupling of z+ and z^ required fi = v. For bo ^ 0, we showed that terms involving [i — v were 
subdominant to the bo terms. Since the same arguments cannot be made for bo =0, we have carried out a dynamic 
renormalisation group treatment in the absence of a mean magnetic field, following the momentum shell approach 
of jl^ . p^ . In our treatment we allow for independent coupling terms XiudxU, and X2bdxb in the u equation and 
Xsdxiub) in the b equation. We assume, as in the scaling argument above, that there is no diagramatic correction to 
the noise strength, and we choose the rescaling so that the noise strength remains fixed. For simplicity we ignore the 
cross-correlation The resulting recursion relations have a stable fixed point at which fi = v. This allows us to 

work with fi — v, and hence to use the decoupled description even when bo — 0. 

We have also carried out a one-loop self-consistent treatment of a slight generalisation of equations (||) and (^) with 
yl„(|fc|) ~ Dik^y^^ , Ab{\k\) ^ Z?2fc '^^^ and fcC(|fc|) — 0. Ignoring noise and vertex renormalisations, we find that 
the fluctuation corrections Sti(fc) and St,(fc) at zero frequency for i> and fi respectively obey: 

Enfc)=/^l ""'P'"^' + (39) 

^' J 27r \p2E"(p)[p2S"(p)H-(fc-p)2S"(fc-p)] ^ ^ 



D2P 



-s+3 



p2Ef'(p) [p^^''ip) + (fc - p)2l]''(fc - p)] 

and 



(40) 



would be of the form (noise-strength/&ofe), while those from the z'^z'^ and z~ z~ would be of order (noise-strength/fc^). The 
former, while in principle dependent on the wavespeed, would be subdominant (since we are now considering z > 1) to the 
latter. 

*This is another point of similarity with 3dMHD: In terms of the Elsasser jl^ variables z*, the nonlinearity in the 3dMHD 
equation for z+ is z~.Vz^, and that in the equation for z~ is z"'".Vz~. The equations in this form have the Galileian invariance 
z* ^ z* -|- Uo, r — > r — Uot for any constant vector Uo, where r is the position coordinate. This guarantees that all the vertices 
in 3dMHD are unaffected by fluctuation-corrections at zero wavenumber. 
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[p^E"{p) + {k - py^"{k - p)] 



(42) 



If we seek a solution of the form E„ — T^k^^ and "Eb = Ti,k^^ we find consistency requires xi — X2 — y/3 = s/3 and 
Tu = r^, i.e., that the scale-dependent kinematic and magnetic viscosities are identical at small k. The singularities of 
the one-loop integral depend strongly upon y. For y > 3, these integrals diverge even when the external wavenumber 
fc 7^ (the "sweeping divergence" For y near 0, on the other hands, one finds S„(fc) ~ ^^(fc) ~ k'^/^. These 

singularity structures are very similar to those seen in a self-consistent treatment of the randomly stirred Navier-Stokes 
equations. 

A technical clarification may be of some interest here. In carrying out the perturbation theory with ho = 0, care 
must be taken to ensure that the condition (6) = is maintained order by order. If ikC{\k\) ^ 0, i.e., there are cross 
correlations in the noise sources in equations (^, and /i 7^ t^, then the perturbation theory generates an apparent non- 
zero < b > and an apparent Alfven wave speed. Including a counterterm to cancel the spurious < b > automatically 
ensures the absence of Alfven waves at bo = 0. At the RG fixed point discussed above, where iJi = v, this issue clearly 
does not arise. 

Our final point concerns the assertion [jllj that the introduction of a nonzero cross-correlation (fufb) in 3(iMHD 
generates, under renormalisation, dissipative terms of the form V^B in the equation for u and V^u in the equation 
for B. If this were true, the renormalised equations of motion would lack the symmetry properties enjoined on them 
by the fact that B and v are a pseudovector and a vector respectively. Unsurprisingly, a straightforward perturbative 
analysis of fluctuation corrections arising from the nonlinearities in (0) rules out any such anomaly. As long as the 



statistics of the noise sources are as in (6a) and (|6c|), so that the transformation properties of v and b are respected, 
no such terms are generated. It is clear that the error in [|ll| arose because their cross-correlation {fu{k)fb{—k)) 
was real and even in wavevector, which is simply not consistent with the nature of the fields in the problem. Such 
a cross-correlation will in fact generate not only the terms mentioned above but terms like dx (ub) in the u equation 
and dx{u^) and dx{b^) in the b equation as well. The resulting system will have only the invariance x —> —x, u —u, 
b — > —6, which has nothing to do with the intrinsic transformation properties of the physical fields. 



V. SUMMARY 

We have shown that a recently proposed ^ Burgers-like one-dimensional model for magnetohydrodynamics 
(IdMHD) is completely equivalent, through a simple transformation, to existing equations in the literature 
We have obtained several new results from these IdMHD equations. The most important of these, which should apply 
to 3dMIID as well, is that the energy spectrum is unaffected by the presence of a mean magnetic field. Apart from 
this, we have presented scaling, rcnormalisation-group and self-consistent analyses of the large-distance, long-time be- 
haviour of correlation functions in the model, and provided arguments for the existence of multiscaling therein when 
the forcing functions have singular small- wavenumber correlations. Lastly, we have corrected an erroneous claim pl| ] 
regarding the effect of cross-correlations between the forcing functions in the velocity and the magnetic field equations. 
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